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We study the formation of Q-balls in the early universe, concentrating on potentials with a cubic or 
quartic attractive interaction. Large Q-balls can form via solitosynthesis, a process of gradual charge 
accretion, provided some primordial charge assymetry and initial "seed" Q-balls exist. We find that 
such seeds are possible in theories in which the attractive interaction is of the form AHtp'tp, with a 
light "Higgs" mass. Condensate formation and fragmentation is only possible for masses in the 
sub-eV range; these Q-balls may survive untill present. 
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I. INTRODUCTION 

Q-balls are lumps of coherent scalar field that can 
be described semi-classically as non-topological solitons. 
They can arise in scalar field theories with a conserved 
global U{1) charge and some kind of attractive interac- 
tion jl]. In a sector of fixed charge, the Q-ball is the 
ground state of the theory. Q-balls generically occur in 
supersymmetric extensions of the standard model Q . In 
these theories, baryon and lepton number play the role 
of conserved charge. 

Q-balls come in two types. Type II Q-balls are associ- 
ated with the flat directions of the potential, which are 
a generic feature of supersymmetric theories. The VEV 
inside the Q-ball depends upon its charge. Formation of 
this type of Q-balls through fragmentation of an Affleck- 
Dine (AD)-like condensate has been studied extensively 
in the literature 1^-^. Type I Q-balls on the other hand 
are characterized by a potential that is minimized at a 
finite VEV, independent of the charge of the Q-ball. We 
have analyzed under which conditions this type of Q-ball 
can be formed in the early universe. In this paper we 
present the results. 

Large Q-balls can form via solitosynthesis, a process 
of gradual charge accretion similar to nucleosynthesis, 
provided some primordial charge assymetry exists [^,0. 
The bottleneck for this process to occur then is the pres- 
ence of initial "seed" Q-balls. Most potentials do not 
allow for small Q-balls which makes solitosynthesis im- 
probable. The exceptions are theories in which the at- 
tractive interaction is provided by a cubic term in the 
Lagrangian of the form AHil;*tp, with a light "Higgs" 
mass. Condensate formation does occur for light fields, 
for masses in the range < eV. If this leads to frag- 
mentation, the thus formed Q-balls can survive evapo- 
ration if their binding energies are large. Q-balls formed 
during a phase transition do not survive evaporation, un- 
less the phase transition takes place at extremely low 
temperatures T < lO^^m^. 

If Q-balls survive untill present they can be part of the 
dark matter of the universe. Recently it was proposed 



that the dark matter could be self-interacting; this would 
overcome various discrepancies between observations and 
predictions based on coUisionless dark matter, such as 
WIMPs and axions Due to their extended nature 
Q-balls have relatively large cross sections, and therefore 
can naturally satisfy the required self-interactions |^ . 

Another cosmologically interesting feature of Q-balls 
is that solitosynthesis in the false vacuum can result in a 
phase transition Accretion of charge proceeds until 
a critical charge is reached, at which point it becomes 
energetically favorable for the Q-ball to expand, filling 
space with the true vacuum phase. 



II. Q-BALLS 

Consider a theory of an interacting scalar field (p that 
carries unit charge under some conserved U{1) charge. 
The potential has a global minimum U{0) — at (p = 0. 
We also require that the potential admits Q-ball solu- 
tions, i.e.. 
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The Q-balls solutions are of the form (p)(x, t) — e'"'(/)(a;). 
The charge and energy of such a configuration is 



Q = Lj d^x (j)- 



and 
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Minimizing the energy for a fixed oj is equivalent to find- 
ing a 3-dimensional bounce for tunneling in the potential 
Uui ■ The bounce solution exists for fiQ < cu < y/U"{0) by 
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virtue of eq. (|l|), and is spherically symmetric The 
equations of motion are: 



r dr 



= 0, 



(5) 



with boundary conditions (f>'{0) — and (f>{oo) ~ 0. 
We will consider scalar potentials of the form 

Ui{<P)= ^mlcj,' - Acj,' + X<l>\ 



(6) 
(7) 



both have 4>o = A/2X and /Xq = — A'^/2\. U2 is 
a typical potential that arises in effective field theories. 
Ui is a non-polynomial potential, as the cube term is 
of the form . It is a typical potential in finite 

temperature theories; this is however not interesting in 
the current context since at high temperatures Q-balls 
evaporate quickly. But it can also arise as an effective 
field theory. Consider for example the potential 



U[{tlj) = m|V> + m]jH*H - A'Hi}*ilj + h.c. 
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where the "Higgs" field H is uncharged under C/(l), 
whereas "0 carries unit charge. Further, we take A' real. 
Now make the field redefinitions 
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then U[ becomes of the Ui form, with (p some particular 
direction in [H, V')-space. We can also calculate /Xg = 
2[/{/(iy9^ cos^ 0) in terms of the U[ parameters. Taking 
rriH — and all quartic couplings equal Ai = A2 = A3 = 
A' to simplify the algebra, this yields 



fields, interacting only gravitationally or through some 
other surpressed interaction with the visible sector. Hid- 
den sectors appear in a variety of models, such as tech- 
nicolor, mirror symmetry, hidden sector SUSY breaking, 
and brane world models. They also arise naturally from 
string theory; in heterotic Eg, x string theory compact- 
ified on a Calabi-Yau manifold, one of the E^s's contains 
the SM whereas the other is some hidden sector. 

We will assume an initial charge asymmetry, i.e., an 
excess of particles over anti-particles. This asymmetry 
may be created through a mechanism similar to those 
invoked to explain the baryon asymmetry in the universe, 
such as the Affleck-Dine mechanism [ |l3[ . 

A. Large Q-balls — thin wall approximation 

For large Q the Q-ball solution can be analyzed using 
a thin wall approximation, which consists of neglecting 
the effect of the surface compared to the bulk. The Q- 
ball may be approximated by a sphere of radius Rq with 
(j) — (j)o inside and zero field value outside. The mass and 
radius of the solition are 



Mq = fiQ, 



and 
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with fi,u! — > fxo for Q — > 00. The soliton is large and its 
cross section is given by the geometrical area 



B. Small Q-balls 



(13) 



3A' 



(10) 



at 6*0 = 7r/4 and Lp^ = 4^73 V2A'. 

Potentials of the form [/{ are present in the scalar sec- 
tor of the MSSM, where the Higgs field couples to spar- 
ticle fields through a cubic interaction . The sparticles 
carry a conserved U{\) charge in the form of baryon or 
lepton number. However, the sparticles and also possi- 
bly formed Q-balls are unstable, as they can decay into 
light fcrmions [0. Stable Q-balls can be obtained in 
a model where the standard model (SM) Higgs field is 
coupled to a charged SM singlet The SM singlet V 
is charged under a hidden sector C/(l)^ global symme- 
try, under which none of the SM particles are charged. 
The Q-balls in this model interact with the SM particles 
only weakly, through the Higgs boson. Another possi- 
bility is that both the H and ip field are hidden sector 



The limit of small charge corresponds to — s- 777^. In 
this limit the solution of the bounce equation (^) is of 
the form |14 
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with a the power of the attractive term in the potential. 
This solution has the right behavior for r — > 00 where 
(/)—!■ and the quadratic term in the effective potential 
dominates, and for uj ~* where the zero of Ui^ is near 
the origin. Using the solution to compute the conserved 
charge (|^), and taking the limit lo —^ m^, one finds a 
finite, non-zero value only for 4 -I- 2D — aD > 0, with D 
the number of spatial dimensions. In Z? = 3 dimensions, 
Ui admits small Q-balls but U2 does not. Therefore, we 
will only consider Ui in the remaining of this section. 

In the limit of large uj, or equivalently very non- 
degenerate minima, one can neglect the quartic terms 
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in C/c^ {(j)) . This is the thick wall approximation ||T^ . The 
approximation is valid for Q-balls with charge Q that 
satisfies, for potential Ui: 



exist for 



Q < U.6m^/y^A, 



Q < 7.3ml/A\ 



(15) 



If above conditions are met one can define an expansion 
parameter 



Q-. 
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with « 4.85. The mass of the soliton is 



Mq = Qm^ 
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The radius of the Q-ball can be parameterized 

rq = 

TO0 



(16) 



(17) 



(18) 



with Pq ~ 0(1). 

The Q-balls described above are classically stable for 
arbitrary small charge Q. However, one expects quan- 
tum fluctuations to become important in this regime. 
Indeed, numerical calculations indicate that this is the 
case, and only configurations with Q ^ 7 are quantum 
mechanically stable . All these calculations are based 
on the potential Ui . This potential is an effective poten- 
tial which is well suited for a semi-classical description 
of large Q-balls. But for small Q-balls the degrees of 
freedom of the underlying potential U[ should be taken 
into account. In this regime a treatment in terms of 
quantum bound states is more appropiate. Solving the 
bound state problem in full generality is not an easy 
task. However in the limit that all quartic interactions 
can be neglected, the theory becomes identical to the 
Wick-Cutkosky model. This model can be solved analyt- 
ically for the case of a massless exchange particle, i.e., 
tuh = 0. The various approaches used in the literature, 
e.g. ladder approximation, Feshback-Villars formulation, 
variational-perturbative calculations |r^-[l9[|, all lead to 
the same result that the bound state spectrum is discrete 
with the n'^ state having an energy (to lowest order in 
a): 



En = 2m^(l 
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(19) 



The above result for the binding energy is derived in 
the limit of a massless boson exchange. No analytic re- 
sults are known for massive scalar exchange. However, 
numerical studies show that bound states still form, pro- 
vided that a is larger than some critical value. We es- 
timate, based on the results in pOl, that bound states 
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That is, the Higgs mass needs to be sufficiently small 
nT-H ^ 1Q~'^{A' /m^Ym^. The energy of the bound state 
is of the same parametric form as for the massless case. 

The other assumption that went into the derivation 
of eq. (|9|) is the absence of quartic couplings. We ex- 
pect this to be a good approximation in the regime where 
quartic interactions are negligible small. The cross sec- 
tion for i/'V'-scattering is a^^^^^ = S\M\^ / IQ-k E^.^. For 
scattering through Higgs exchange, governed by the cu- 
bic interaction, this gives at tree level 
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Here E = maxjT, m/f }. At low temperatures T < A', 
which are the temperatures of interest, the cross section 
quickly approaches the unitarity bound and higher or- 
der diagrams cannot be neglected. In this regime we 
will approximate the cross section by cr ^ ""^^ with 

— 27: /m^ the Compton wavelength. Scattering 
through the quartic point interaction has an amplitude 

— A'. And thus the requirement that the repulsive 
quartic interactions are negligible small o-quartic ^ ^cubic^ 
are fuUfiUed for all quartic couplings A' ^ 1. It may be 
that also for non-perturbative values of the quartic cou- 
plings bound states persist; but this certainly cannot be 
analysed perturbatively. As it seems unnatural to have 
quartic couplings larger than one, we will ignore this pos- 
sibility. 

On the other hand, the quartic couplings cannot be ar- 
bitrary small or else no Q-ball solution exists: for the case 
of zero Higgs mass and all quartic couplings equal, /Iq > 
translates into A' > A'^/3m^, as follows from eq. (|l0|). ^ 
Non-zero, but small Higgs mass mu < lO^^m^ does not 
alter this result noticeably. The quartic couplings do not 
have to be all equal, but at least one of them has to 
be 0{A' /m^y' . For A' = m^, Q-ball solutions exist for 
example for (Ai,A2,A3) = (0.4, 0.4, 0.4), (1, 0.01, 0.01) 
and (0.05, 0.8, 0.3). 

Both the quantum bound states discussed above and 
Q-balls describe the same objects — stable bound states 
with a fixed global charge — but in a different language. 



^Condition /ig > corresponds to the requirement that <j) = 
is the global minimum of the potential. Q-ball solutions do 
exist for = a local minimum. In the potentials Ui and U2 
this possibility is not realized, since at low temperatures the 
field will end up in the true vacuum. {Ui: at the temperature 
T that the minimum at 7^ becomes global the energy 
barrier is ~ lO^^AT*. U2: at high temperature m?{T) < 0.) 
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In both descriptions the existence and stabihty of the 
bound state reUes on the trihnear couphng and the con- 
served global charge (that is, conserved particle num- 
ber). For large bound states quantum fluctuations can 
be neglected, and a semi-classical description as a Q-ball 
becomes a good approximation. The trilinear coupling 
makes it possible for the energy of a bound state with a 
fixed charge to be less than a collection of free particles 
with the same charge. In the limit of small particle num- 
ber (global charge) , it becomes necessary to treat the full 
quantum problem because the semiclassical approxima- 
tion breaks down. The trilinear term can be viewed as 
an attractive interaction between the 0-particles, which 
makes it possible for bound states to form. The lowest 
level bound state is the stable ground state, as charge 
conservation forbids it to lower its energy through anni- 
hilation of (j) particles. 

It is tempting to compare the ground state result 
(n = 1) of eq. ( |l9[ ) with the q = 2 result obtained in the 
thick wall approximation (|l7|): both mass formulas give 
the same parametric dependence. However, in the over- 
lapping regime both approximations are taken beyond 
their domain of validity: for equal masses mu = 
bound states can only form for large a, and for q = 2 a 
semi-classical treatment breaks down. Of course both ap- 
proximations are similar in that they neglect the quartic 
interactions. 

In conclusion, the potential U[ admits stable, two- 
particle bound states at low temperatures (below the 
binding energy), provided the Higgs mass is sufficiently 
light, and the quartic repulsive interactions small. We 
repeat that our assumption here is that non-zero quar- 
tic couplings do not destabilize the bound state provided 
< <^{^^t4,4,)i this should be checked by an ex- 
plicit calculation. For the potential to have a global min- 
imum at = 0, or equivalently for Q-ball solutions to 
exist into which the bound states can grow, the couplings 
cannot be too small. 



A' < 1 



A' > 



A'' 



repulsive forces small 
small Higgs mass (^22) 
Q— balls exist 



A possible set of parameters is A' ~ 0.5, A' ~ and 
iTiH ^ lO^^m^. The binding energy for the bound state 
is then B2 = ~ 5 x lO^^m^, and /io ~ 0.6?7i^. 

We will further assume that similary bound states of 
more than two particles can exist, and that they have 
energies 



Mr 



a 



1 A'^ 



(23) 



with /q some unknown factor depending on the charge 
Q, the mass of the exchange particle and the strength of 
the quartic interactions. 



The binding energy of a Q-ball is Bq = Qm^ — Mq. 
Two-particle bound states are only stable at tempera- 
tures below the binding energy T < B2 ~ a^/S. From 
then on they can grow by capturing charged particles. A 
non-relativistic particle with kinetic energy ^ T has 
energy T -\- Bq inside the Q-ball/bound state. For it to 
be captured it has to lose an amount larger than T in 
the collision. Assuming isotropy, on average a particle 
will lose half of its energy. Therefore, for temperatures 
T < Bq a considerable amount of the particles scattering 
with the Q-ball will be captured. We approximate the 
absorption cross section o-abs(Q) for a Q-ball with charge 
Q by the scattering cross section: (Tabs(Q) ^ ""^q- 

III. SOLITOSYNTHESIS 

In thermal equilibrium, the production of large Q-balls 
through gradual charge acretion is very efhcient. This 
process is dubbed solitosynthesis for its similarity with 
nucleosynthesis. It requires an initial charge asymmetry 
not unlike the baryon asymmetry of the universe. Freeze 
out of any of the reactions involved will put a halt to 
solitosynthesis. 

In this section we will describe the thermodynamics of 
Q-balls in terms of a gas of non-relativistic '0 particles in 
thermal equilibrium. The '0 particles can bind together 
through the exchange of a light scalar particle, as given 
by the cubic interaction in U[. For large Q-balls a semi- 
classical description in terms oiUi suffices, and can be 
replaced by 4> in s\\ the formulas. 



A. Q-balls in thermal equilibrium 

At non-relativistic temperatures T < m^, the number 
densities of a distribution of Q-balls and free ip particles 
in thermal equilibrium are governed by the Boltzmann 
distribution: 



and 



nQ{T) ^gQ 



MqT 

1-K 
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^{pq-Mq)/T 



(24) 



(25) 



Here is the internal partition function of the Q-ball, 
and 5^ = 2, the number of degrees of freedom of a com- 
plex field. Solitosynthesis is only possible if capture rates 
are large compared to the expansion rate of the universe, 
otherwise the densities are frozen. If so, the gas of ?/; 
particles and Q-balls is in chemical equilibrium, and the 
accretion and absorption reactions 



(Q) + v^(g + i) 



(26) 
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enforce a relation between the various chemical poten- 
tials: fiQ = Q^ip. This allows one to express the Q-ball 
number density in terms of the i/i-number density 



nQ{T) 



9Q_ 

Q 



Mq 
my, 



3/2 
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3(Q-l)/2 



.Bq/T 



(27) 



with Bq = Qm^, — Mq > the binding energy of a 
Q-ball. Similar equations can be written down for the 
number densities of anti-i/^'s and ant i-Q- balls. 

We will assume a primordial asymmetry of ip^s over 
^*'s, T] = (n^ — n^')/n^, where = 2.4T^/7r^ is the 
photon number density. Initially one has both "0 and ip* 
particles. For large asymmetry the anti-particle density 
can be neglected. Also, if the Higgs mass is light then pair 
annihilation occurs, and at non-relativistic temperature 
anti-particles deplete rapidly. The annihilation reaction 
enforces — —^^p-, which in the non-relativistic limit 
leads to 



(28) 



For temperature T ^ the chemical potential /i ; 
otherwise the Boltzmann suppression exp[(/i — m)/T] is 
tremendous and the charge conservation equation 



(29) 



can never be satisfied. Annihilation is efficient untill the 
density of anti-particles is negligible small. The number 
density of stable V'-particles is then 



2.5 X IQ-^VL^h^ 



GeV 



(30) 



The local density can be higher if clustering occurs. This 
is not to be expected until the matter dominated era, 
T < Teq « 5.5(f}o/i^)"^C"^eV. 

The photon temperature may in general be different 
from the temperature of the Q-ball system. Particle 
species that decouple from the heat bath when they are 
highly relativistic maintain an equilibrium distribution 
with temperature T oc . The photon temperature red 
shifts as cx g^s^^R~^, and thus the difference in tem- 
peratures is given by a factor ( = {g^^siTo) — g*s{T)Y^'^ , 
with Td the temperature at which the Q-ball system 
decouples. When the ^/i-particles only interact gravita- 
tionally C ~ 10, whereas it can be much lower for more 
general interactions < — > X, where X are light par- 
ticles that do not carry the same Uq{1) charge as the 
■0-particles. We parameterize = (T, with C ~ 1 — 10. 

The Q-ball densities can start growing when the ex- 
ponent in equation ( p7| ) dominates over the potentially 
small factor in front. Since Bq grows with Q, forma- 
tion of large Q-balls is favored. The evolution of a single 
Q-ball is given by the absorption and evaporation rates 




FIG. 1. Accretion time Tg plotted as a function of cC^'^rj 
for various values of IqIvti^. For large values of cC^^t) — to 
the left of where the lines stop — accretion dominates over 
evaporation at all (non-relativistic) temperatures 



of -0 particles by a Q-ball of charge Q [^. These can 
be found using detailed balance arguments. In chemical 
equilibrium we have for the process in eq. (Eq) 



"V'^'/' f^abs(Q) = "-Q+1 ?'ovap((3 + 1), (31) 



and 



^ ^abs(Q) ^cvap(Q) 



n^v^ [0'abs((3) 



nQ-\ 
nQ 



fTabs(g-l)]. (32) 



The Q-ball starts growing when Tabs (Q) > f'ovap((3)- This 
happens for T ^ Tg with 



T 



(2^)-ln(cC3r;)' 



■In 



and 



Iq = + Mq^i - Mq. 



(33) 



(34) 



Here c = {Q/Q — lY^^^^^ gg/gQ-i, which goes to one 
for large Q, and c ~ 10 in the limit Q ^ 2. Iq is the 
binding energy with which a single -^-particle is bound 
to the Q-ball. For very small Q-balls /g/m^, = 
Accretion of the smallest (QmiiO-ball starts when T ^ 
-^Qmin < lO^'^TO^ for A' < 1. For large Q-balls Iq = 
TT^i! ~ ^^o■ Figure 1 shows Tg as a function of cC^t] for 
various values of Iq. For large equation (|3^) has no 
solution; here absorbtion dominates at non-relativistic 
temperatures. 

For large H mass and early freeze out the charge asym- 
metry may be small, as then both ijj and ■i/;*-densities are 
large at freeze out while annihilation is negligible. Both 
particle and anti-particle number are conserved, and one 
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can in principal have growing Q- and anti-Q-balls at the 
same time. In this case however, the formation of small 
seed Q-balls, which are necessary to start the fusion pro- 
cess, appears to be a major obstacle. 



B. Freeze out 

For solitosynthesis to work Tg must be higher than 
Tfo, the temperature at which the absorption reactions 
eq. ( |2^ ) freeze out. This occurs when the reaction rate 
for accretion becomes smaller than the expansion rate of 
the universe: 



T[{Q)+i,^{Q + \)]<H. 



(35) 



The Hubble constant during the radiation dominated era 
is H = 1.7gli^T^/Mpr, it is if = 1.7 gli''T^^^Ta,i^ /Mpi 
in the matter dominated era. The effective degrees of 
freedom are g^s < 10 for T-^ < lOMeV and g^s 100 
for O.lGeV < < lO^GeV. The accretion rate is 
r — n^t;^,f7abs((3)- Neglecting the charge density sub- 
siding in Q-balls, then for stable tp particles is given 
by eq. (|3G|). At non-relativistic temperatures = 
(r/2^m^)V2. 

We are interested in temperatures T < Tg] then the 
cross section is (Tabs ^ '^^q ^-^^ freeze out occurs for 
temperatures, for Tpo > « 5.5{floh'^)~^(~^eV: 



T, 



FO <; 
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And for Tpo < TLq 



< 



lO'^-^ / Y 



3/2 / 0.3 \ ( qV^ 
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(36) 
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For A' ~ m^, C 10 and < GeV freeze out of the ac- 
cretion reactions for the smallest Q-balls {Q — 2) occurs 
after the accretion phase, Tpo ^ Tg. In this case soli- 
tosynthes can start at photon temperature Tj = QTg ~ 
IQ-^rn^ < lO^'eV. Note that Iq cx [A' /m^)'^ decreases 
rapidly for smaller cubic couplings and Tpo < Tg can 
only be satisfied for increasingly low -(/'-mass. In the mat- 
ter dominated era the reaction rate, and thus the freeze 
out temperature, can be increased through clustering. 
For an overdensity of ^ 10^ in galaxies, we find that for 
A' ~ O.Ito^ solitosynthesis occurs for m,0 < GeV, start- 
ing at temperatures T^ ^ eV. For smaller values A' < 0.1 
or ^ MeV, solitosynthesis occurs in the future, at 
temperatures smaller thant the present day temperature 
T-y <Tq^ 2.35 X 10-4 eV. 

With the accumulation of charge in Q-balls the number 
density of V' particle decreases and the system freezes out. 



Since the accretion is such an explosive process, this will 
generelly not happen until almost all charge resides in Q- 
balls. More quantitatively, when the ■i/'-density decreases 
to 10% of its original value, Tpo decreases only by a factor 
2Q-2/3 rpj^g back reaction is only important when Tg « 
Tpo, and it shuts off the growth of Q-balls immediately; 
in all other cases most charge will end up in Q-balls. 

The accretion rate of a single Q-ball is limited by the 
diffusion rate. However, diffusion of charge is only im- 
portant when / ^ Rq, with I ~ F^^ the mean free path. 
The radius Rq of a Q-ball becomes equal to the mean 
free path for a large charge: 



QdiffCT) - 10 



87 



m,p/T \ 
J 



7/2 



1 0.3 



(3 nh^ GeV 



(38) 



For Q > Qdiff diffusion is important. The total amount 
of charge inside a Hubble volume is Qtotai = n-^H^^ 



Qtotal(r) - 10' 



63 / m^/T 102 



107 



Cg 



1/2 



GeV 



0.3 



(39) 



For small masses Qdiff may be lower than the total charge 
inside a Hubble volume; in this case Qdiff will be an up- 
per limit on the charge of the Q-balls formed during soli- 
tosynthesis. 



IV. SEEDS 

Solythosynthesis is a very effecient way to form large 
Q-balls, provided there are some initial seed Q-balls at 
temperatures above freeze out. These seeds may be rem- 
nants of an earlier epoch, formed during a phase tran- 
sition or via the decay of a Bose-Einstein condensate. 
Another option is that small stable Q-balls can form in 
the gas of (/)-particles. 



A. Formation of small Q-balls 



As discussed in section HE , small Q-ball solutions are 
only stable for potentials with a cubic interaction. Two- 
particle bound states can form through scalar exchange, 
provided the mass of the exchange boson is sufficiently 
small and the quartic interactions can be neglected. In 
this case seed Q-balls can be formed copiously and soli- 
tosyntesis can start. If the mass of the scalar mass is 
of the same order as the mass of the charged particles 
two-particle bound states do not form, but it may still 
be that small Q-balls with charge Qmin > 2 are stable. 
Numerical calculations indicate that in the thick wall ap- 
proximation (which has mn = ni^) Q-balls are quantum 
mechanically stable for Qmin ^ 7 p^. 
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For Q > 2, Q-ball formation is surpressed compared 
to the two-particle bound state, by the requirement that 
Q charges should be similtaneously in a volume of radius 
~ Rq. Define P{q) to be the probability to find a charge 
Q in the volume of a Q-ball, Vq ~ R^- The mean charge 
in Vqisq — n^Vq, whereas the variance is — ((Ag)^) — 
T{dq/dfi)T.v = 9- Since 



9« l.OCV 



T 



< 1 



(40) 



a discrete distribution is needed, the Poisson distribution: 



q 



(41) 



The density of lumps with charge Q in a volume Vq is 
Uq — P{q)/Vq. The reaction rate for the bound state 
is Fg"'^ ^ a^^'^nq so that the chance that in a Hubble 
time a "Q-lump" forms a bound state is ~ nqO^^'^H^^ . 
Multiplying this with the total number of Q-lumps in a 
Hubble volume gives the number of Q-ball seeds Nq ^ 
nla^q'"^H-^. Taking Rq - this yields 



iV,^^«(ab„d™l) — -r^ • (42) 



Assuming Cg"'^ ^ this gives an upper bound on Nq. 
Only for small q = 2, 3,4 or so Nq is larger than unity, 
and there is seed forming. 



B. Primordial seeds 

The seed Q-balls may also be Q-balls formed at an 
earlier epoch. For this to be possible the initial Q-balls 
should be large enough to survive the period of evapora- 
tion. The evaporation rate is given by the detailed bal- 
ance equation (^). Ignoring absorbtion, which is sub- 
dominant for T < Tg (note that the evaporation rate 
decreases exponentially with temperature) , one gets that 
the smallest Q-ball to survive has charge 



Qs - 10" 



13^ /GeV 




dT^^ 

Iq/t 



(43) 



with Ti the temperature at formation. For Tpo, Tg < Ti, 
the integral can be approximated by ~ exp(— /q /Ti) + 
Iq /Ti). Only for masses ^ eV is Qs smaller than the 
total number of particles available in a Hubble volume 
at Ti ~ m^, eq. (p9|), and is there a change for very 
large primordial Q-balls to survive the period of thermal 
evaporation. 

Another possibility is that formation happens at the 
onset or during the acretion phase: ^ Tg. For large 
binding energy Iq — > (which is possible for large Q- 
balls) and large r/, accretion dominates over evaporation 
at non-relativistic temperatures, see figure 1, 



Primordial Q-balls may also form during a first order 
phase transition from the false "Q-ball vacuum" to 
the true vacuum. At the Ginzburg temperature ther- 
mal transitions between regions of false and true vacuum 
freeze out; any region of false vacuum with a charge larger 
than the minimum charge of a stable Q-ball surviving 
below this temperature will become a Q-ball. The po- 
tentials under considerations do not exhibit the required 
first order phase transition (see footnote 1). One could 
add additional terms to the potential to get a phase tran- 
sition. However, the survival of regions of false vacuum is 
exponentially surpressed with size, and correspondingly 
Q-ball formation is exponentially surpressed with charge. 
If formed, the Q-balls are expected to be small Q ~ Qmin- 
Unless there is a mechanism to delay the phase transi- 
tion to very low temperatures T ^ 10~^m^, these Q-balls 
quickly evaporate and are cosmologically unimportant. 

Formation of primordial Q-balls through fragmenta- 
tion of a condensate 101 is studied in the next section. 



V. BOSE-EINSTEIN CONDENSATION 

We will now study whether there will be condensation. 
A condensate that is unstable under fluctuation can frag- 
ment into possibly large Q-balls. We will consider the 
effective potentials Ui and C/2- In this section — 1, 
i.e., all quantities are expressed in units of mass. 

We will assume that the number density of anti- 
particles can be neglected and p « n^. The state of the 
system is given by the minimum of the effective potential 
for a fixed charge Q: 



V{q,^)^V{n,<j>)+fip, 



(44) 



with V"(/i, (/)) the effective potential for a fixed chemical 
potential. In this section denotes the classical back- 
ground field, and (po its value at the minimum of V{q, (j)). 
A non-zero value of 0o signals the existence of a conden- 
sate. At finite temperature the freqency lu of the Q-ball 
can be identified with the chemical potential ^ . The 
charge density can be solved from 



d/i 







P = P{p)- 



(45) 



Eliminating /i in cq. (^4|) then gives the effective poten- 
tial in a fixed charge section. A stable configuration lies 
at the minimum of V{q, (j)). 

To analyze the stability of the condensate one can con- 
sider fluctuations in the homegeneous background. From 
the dispersion relation it follows that fluctuations are am- 
plified for wavelengths smaller than fcmax P,p2[: 



u"{4>o)- 



(46) 
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For — (j)QU"{(f)o) < the above equation does not have 
a physical solution and the condensate is stable. 
We parametrize the charge density is p = rjUy 



3 X IQ-^C^ 



0.3 



GeV \ 
J 



rp6 



ri'T^. (47) 



A. Non-Relativistic Limit 

At zero temperature the charge density p oc T'^ is zero, 
and there is no condensate. At non-zero temperature 
condensate formation will occur if the charge is larger 
than the number of excited states. 

In the non-relativistic limit the finite-temperature cor- 
rections to the potential are small, and as a first ap- 
proximation we can use the zero temperature result 
V{n,(j)) = U{(t)) - l/2pV^ with U{(t)) the classical po- 
tential, eq. (||, ^. Equation ( ^ ) gives p — fic/)'^. The 
minimum of V{q, <j)) is at 



pI 



4>l 

6^ 



^3Acf>l 



for Ui 
for C/2- 



(48) 



At low temperatures p ^ 1 and p w 0q <C 1. In this 
limit a possible condensate is unstable against decay for 
values (/)o < 14 for Ui and 0o < for U2, as follows 
from eq. (|4^). 

To see whether a condensate actually forms one has 
to compute the density of thermal states. At low tem- 
peratures the cubic and quartic terms in the potential 
become negligible small, and the theory approaches the 
free theory. In this limit the number of thermal states is 



C(3/2) 
(27r)3/2 



T 



3/2 



(49) 



Since p oc T^, at low temperature all charge will be in 
the excited states and the condensate is empty. The only 
chance to have a filled condensate is for T ^ 1 and rj' 
large, so that p > or r^'T^/^ > C(3/2)/(27r)3/2 _ q 
Note however that in the limit T — > 1 the non-relativistic 
approximation breaks down, whereas in the limit rj'T^ — > 
1 the free field approximation breaks down. 



B. Relativistic limit 

We will first consider the potential Ui. The effective 
potential for fixed chemical potential to highest order in 
T is 



y(p,0) = 1/2(1 + ATVS-p 
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c(T) + 0(T), 



(50) 



with c(T) some temperature dependent constant which 
we will drop. From this it follows that 



p = pcj) + pT 



(51) 



The first term in the above equation is the charge in 
the condensate, the second term represents the charge in 
excited states. The charge fraction in the condensate is 
4>l/{(f>l + T^), which is small for 0o < T. The efi'ective 
potential for fixed charge density in the relativistic limit 
becomes 



A, 



A(tf + Ac/)'' 



3p2 



2(3(/)2 + T2) ■ 
(52) 



Consider the case 0o ^ T] then the potential is mini- 
mized at ^0 = (and thus the approximation is consis- 
tent), provided 



,2 A 1 
^ 27 ■ 



(53) 



This can also be seen from the second derivative 
F"(g,0) = 1 + A/3T2 - 97?'^T2, which becomes nega- 
tive for large rj' . Thus if condition ( p3| ) is obeyed there 
is no condensate. For finetuned values of ^2/;^ ^ second 
minimum of the potential may develop, but since in the 
limit of large temperature the only minimum is at (/)o = 
the field will not end up there. 

Consider then the potentially more interesting case 
that 1]' is large, and condition ( |53| ) is not satisfied. Then 
V"{Q) < and the potential is minimized at non-zero 
field value. Minimization in the limit T ^ ipo as well as 
in the limit T <C (/lo does not yield a consistent solution. 
It follows that the minimum is at i^o ~ This is con- 
firmed by numerical calculations. The charge density in 
the condensate is comparable to that in excited states. 
The condensate is unstable for /Cj^ax > O7 6q. (HQ), with 



2y6 



1 



A(r2 -|-360g) - 1. (54) 



At large temperatures (j>o oc A""' and the second term 
in the above equation dominates, as can be verified nu- 
merically. The condensate is stable for large T. The 
condensate becomes unstable in the limit T 1. As this 
is also the limit in which the high temperature expansion 
breaks down, it is unclear whether the condensate really 
fragments. 

The analysis for potential U2 is similar. At high tem- 
perature the effective potential becomes 



V{p,c^) = lil-lAT'-p')f 



-c{T) + 0{T). 



(55) 
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For this case, the equivalent of eq. (p3[) is 



V 



A_ J_ 
12 97^' 



(56) 



At large temperature a stable condensate will form with 
(j)o \/ A/2 + 6r]'^/2VX for small asymmetry or (f>o 
T for large asymmetry rj' > 1/9. The condensate only 
survives in the limit T — > 1 for large t]' . The condensate 
may be unstable in this limit. 

To conclude this section, at non-relativistic temper- 
atures there is no condensate and all charge resides in 
excited states. At temperatures T > consensation 
occurs for large assymetries t]' > 1/9, corresponding to 
masses m,^ ^ eV. The condensate becomes unstable in 
the limit T — s- and fragments into Q-balls. Cau- 
tion should be taken, as he high temperature expansion 
breaks down in this limit. If the binding energy of the 
Q-balls is sufficiently large /2\ > 10~^ the period of 
evaporation is absent, see figure 1, and these Q-balls sur- 
vive. 



VI. CONCLUSIONS 

To summarize, solitosynthesis is a very effecient way 
to form large Q-balls provided some primordial charge 
asymmetry and initial seed Q-balls exist. Most theories 
do not allow small stable Q-ball or bound state solutions, 
and solitosynthesis does not start. The exception are the- 
ories in which the attractive interaction is provided by a 
cubic term of the form AH^j*ip. Bound states can form 
if the Higgs mass is light [mH/m^ < lO^'^A^/m^). No 
bound state calculations have been done in the presence 
of quartic coupling. We assume that for quartic interac- 
tions that are small compared to interactions governed 
by the cubic term bound states persist (A ^ 1). We note 
that if this assumption turns out to be too optimistic, 
and stable bound states require smaller quartic coulings, 
then small bound states and large Q-balls become mutu- 
ally exclusive. This is because for the potential to admit 
Q-ball solutions the quartic coupling cannot be too small 
(A > /rn^j). SuccesfuU solitosynthesis will occur if 
the accretion phase happens before the system falls out 
of equilibrium. All these conditions together limit the 
paramater space severely. 

For solitosynthesis to have happened in the early uni- 
verse one needs A = 0.1 — Im^, at least one of the 
quartic couplings A ~ 1, mn ^ lO^^m^, and MeV < 
^ GeV. This rules out models in which the H field 
is the standard model Higgs field, such as the MSSM 
and the model studied in p^ . The temperature at 
which Q-balls start growing decreases very rapidly with 
A: Tg/m^ cx (A/rn^)**. For smaller values of the masses 
or of the cubic coupling than given above, solitosynthesis 
may still happen in the future. 



Q-balls that survive untill present can be part of the 
the dark matter in the Universe. For them to play a role 
during structure formation they must have been formed 
before the universe became matter dominated, that is at 
temperatures > Toq — 5.5(r2o/i^)^^eV. This is only 
possible for A' ^ 1 and ~ GeV. Whether the Q- 
balls can fuUfill the required cross section to mass ratio 
to overcome the problems with cold dark matter as pro- 
posed in 1^ remains another question. More (numerical) 
studies are needed to determine if solitosynthesis results 
in a few Q-balls with a very large charge, or in a large 
number of Q-balls with lesser charge. 

Condensate formation is only possible for large asym- 
metries, or equivalently ^ eV. Symmetries of 
the order one can be generated through the Affleck- 
Dine Mechanism ||l^. Early decoupling increases the 
number of charged particles by a factor C,^ with C, = 
{g*s{TD) — g^:s{T)Y/^, which favors condensation. The 
condensate becomes unstable against fluctuations in the 
limit T m^, i.e., the limit in which all the used ap- 
proximations break down. Evidently, better approxima- 
tions are needed to settle the matter. Q-balls formed 
through a possible fragmentation of the condensate sur- 
vive untill present if accretion dominates over evapora- 
tion at non-relativistic temperatures. This is possible for 
Q-balls with a large binding energy, Iq = — (m^ — 
AV2A)i/2 > 10-2^^. 

The potentials studied do not allow for a first order 
phase transition from the false "Q-ball vacuum" to the 
true vacuum. One could try and add terms to the po- 
tential so that such a phase transition occurs. However, 
the Q-balls that may form during the phase transition 
are small and will evaporate quickly. 

Solitosynthesis can lead to a phase transition from the 
false to true vacuum. This will not happen for the po- 
tentials studied in this paper, as for these the field will 
always end up in the true vacuum. 

The author thanks Alexander Kusenko for very helpful 
discussions. This work was supported in part by the US 
Department of Energy grant DE-FG03-91ER40662, Task 
C, as well as by a Faculty Grant from UCLA Council on 
Research. 
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